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c — c correlations are calculated in the fct-factorization approach. Different unin- 
tegrated gluon distributions (uGDF) from the literature are used. The results are 
compared with recent results of the FOCUS collaboration. The recently developed 
CCFM uGDF gives a good description of the data. Predictions for the HERA 
energies are presented. 

1 Introduction 

In recent years a lot of activity was devoted to the description of the photon-proton 
total cross section in terms of the unintegrated gluon distribution functions (uGDF) 
(see e.g. |H2| and references therein). In some of the analyses also inclusive charm 
quark (or meson) were considered 0. Although the formalism of uGDF is well 
suited for studying more exclusive observables, only very few selected cases were 
considered in the literature. A special example is azimuthal jet-jet correlations 
in photon-proton scattering 0|. Recently the FOCUS collaboration at Fermilab 
provided precise data for c — c correlations [Sj- We analyze the c — c correlations 
in terms of uGDF available in the literature. While the total cross section depends 
on small values of x, the high-pt jets and/or heavy quark production test gluon 
distributions at somewhat larger x. Only some approches are applicable in this 
region. In particular, we test results obtained with recently developed CCFM 
unintegrated parton distributions. This presentation is based on where more 
details can be found. 

2 Charm-anticharm correlations 

The total cross section for quark-antiquark production in the reaction 7 + p — > 
Q + Q + X can be written as 0] 

a^QQ(W) =jd<t>j dpi, J dpi, J dz /g( y 2) • a{W,Pi..uP2,u z) . (1) 

In the formula above f g (x,K ) is the unintegrated gluon distribution. The gluon 
transverse momentum is related to the quark/antiquark transverse momenta pij 
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and p2t as: 



k 2 = v\,t+v\.t + 2piAP2,tCOSC, 



(2) 



Above we have introduced: 
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The unintegrated gluon distribution f g is evaluted at 
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The basic ingredient of our approach are unintegrated gluon distributions. Dif- 
ferent models of uGDF have been proposed in the literature (see for instance |1I2| ). 
The main effort has been concentrated on the small-x region. While the total cross 
section is the genuine small-x phenomenon (x < 10 -3 ), the production of charm 
and bottom quarks samples rather the intermediate-x region (x~ 10~ 2 - 10 -1 ) even 
at the largest available energies at HERA. It is not obvious if the methods used 
are appropriate for the intermediate values of x. In the present approach we shall 
present results for a few selected gluon distributions from the literature. For il- 
lustration we shall consider a simple BFKL JT], saturation model used to study 
HERA photon-proton total cross sections ^2] (GBW), saturation model used re- 
cently to calculate particle production in hadron-hadron collisions JHj (KL). These 
three model approaches are expected to be valid for low values of x. At somewhat 
larger values of x all these models are expected to break. As in Ref.^Hl we shall 
try to extend the applicability of the model by multiplying the model distributions 
by a phenomenological factor (1 — x) n . 

The azimuthal correlation functions w((j>) defined as: 



w{4>) = / dp l t / dp; 



At 



dz 



fg(Xg, K ) 



v(W,pl,t,P2,t,z) . 



(G) 



and normalized to unity for two energies of W = 18.4 GeV (FOCUS) and W — 
200 GeV (HERA) are shown in FigHJ The GBW-glue (thin dashed) gives too 
strong back-to-back correlations for the lower energy. Another saturation model 
(KL, Q|) provides more angular decorrelation, in better agreement with the ex- 
perimental data. The BFKL-glue (dash-dotted) provides very good description of 
the data. The same is true for the CCFM-glue (thick solid) and resummation-glue 
(thin solid). The latter two models are more adequate for the lower energy. The 
renormalized azimuthal correlation function for BFKL, GBW and KL models are 
almost independent of the power n in extrapolating to larger values of x g . In cal- 
culating the cross section with the CCFM uGDF for simplicity we have fixed the 
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scale for u 2 = 4 m^. In the present calculations we have used exponential form 
factor with b e — 0.5 GeV -1 (see [Rj\). For comparison in panel (b) we present 
predictions for W = 200 GeV. Except of the GBW model, there is only a small 
increase of decorrelation when going from the lower fixed-order energy region to 
the higher collider-energy region. 




Figure 1. Azimuthal correlations between c and c. The theoretical results are compared to the 
recent results from |S| (fully reconstructed pairs). 



In the LO collinear approach the transverse momenta of two jets add up to 
zero. It is not the case for our LO k t -factorization approach. In Fig|3|we present 
normalized to unity distribution in p^_, where p+ = Pi + P2- Due to momentum 
conservation, in our approach the sum of transverse momenta is directly equal to 
transverse momentum of the gluon {p + — k) . This means that the distribution in p^_ 
directly probes the transverse momentum distribution of gluons. The CCFM gluon 
distribution gives the best description of the FOCUS data We expect that this 
approach is suitable for x > 0.01. Although the other models give also reasonable 
description of the FOCUS data, one should remember that their application for the 
low energy data is somewhat unsure. 

Not only azimuthal correlations are interesting. The formula Q can be rewrit- 
ten in the form 

a^QQ {W ) = J dpi J dp 2 2tt w(pt t , P l t ; W) , (7) 

where the two-dimensional correlation function 

w{pl tl p\ t - W)=Jdcb dz /g( y 2) • *(W,Pi,t,p 2 ,t, z) . (8) 

In Fig|2we present two examples of w(p\ t ,p\ t ) at W=18.4 GeV for different uGDF. 
As in Ref.^Hl f° r GBW we have used the power n=7 in the extrapolating formula. 
While normalization is slightly dependent on the value of the power, the shape of 
the two-dimensional map is almost the same. The maps for different uGDF differ 
in details. The distribution for the GBW gluon distribution is concentrated along 
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Figure 2. The two-dimensional maps w(p\ t ,p| t ) for W = 18.4 GeV: a) GBW and b) CCFM. 

the diagonal p\ t = p% t and in this respect resembles the familiar collinear LO 
result. The CCFM distribution has a sizeable strength at the phase-space borders 
for p\ t or p\ t ~ 0. Experimental studies of such maps could open a possibility 
to test models of uGDF in a more detailed differential way. In principle, such 
studies will be possible with HERA II runs at DESY. 




Figure 3. distribution of c — c. The theoretical results are compared to the recent results from 
5 (fully reconstruced pairs). The notation is the same as in FigUl 



3 Summary 

We have shown that analysis of kinematical correlations of charm quarks opens a 
new possibility to verify models of uGDF. The recently measured data of the FO- 
CUS collaboration at Fermilab allows to study the unintegrated gluon distribution 
in the intermediate-x region. The recently developed unintegrated gluon (parton) 
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distributions which fulfil the CCFM equation describe the data fairly well. Many 
models of uGDF from the literature are constructed rather for small values of x 
and its application in the region of somewhat larger x (x > 0.05) is questionable. 
It can be expected that the correlation data from the HERA II runs will give a 
new possibility to verify the different models of unintegrated gluon distributions in 
more detailed way. 
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